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Abstract 

The angular power spectrum Ci is extracted from the 4 year COBE 
DMR data with a 20° galactic cut, using the narrowest window func- 
tions possible. The average power in eight multipole bands is also 
computed, and plotted together with a compilation of power spectrum 
measurements from other experiments. The COBE results are found 
to be consistent with the n — 1 power spectrum with the normaliza- 
tion Q„ns,ps = 18 /iK reported by the COBE DMR team. Certain 
non-standard cosmologies, such as "small universe" models with non- 
trivial spatial topology, predict power spectra which are not smooth 
functions. Rather, they contain bumps and wiggles that may not have 
been detected by other data analysis techniques such as the Hauser- 
Peebles method (Wright et al. 1996), the band-power method (Hin- 
shaw et al. 1996) and orthogonalized spherical harmonics (Gorski et 
al. 1996), since these methods all give broader window functions that 
can smear such features out beyond recognition. On the large angular 
scales probed by COBE, the Universe thus appears to be kind to us, 
presenting a power spectrum that is a simple smooth function. 
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1 INTRODUCTION 



Since their discovery (Smoot et al. 1992), fluctuations in the Cosmic Mi- 
crowave Background Radiation (CMB) have emerged as possibly one of the 
most promising ways of measuring key cosmological parameters such as the 
Hubble constant h, the density parameter 0, the baryon fraction f^, the 
cosmological constant A, etc. As the angular power spectrum of the CMB, 
usually denoted Cg, depends on about a dozen cosmological parameters (see 
e.g. White et al. 1994, Bond 1995a, Hu 1995, Steinhardt 1995 or Tegmark 
1996b for recent reviews), accurate determination of this power spectrum 
by future experiments could be used to measure all these parameters simul- 
taneously, with accuracies of a few percent or better (Jungman et al. 1996). 
Although we are still years away from attaining this goal, which would re- 
quire mapping the CMB at high resolution over a large fraction of the sky, 
our knowledge of the power spectrum has been growing steadily over the last 
few years. The first year of COBE/DMR data (Smoot et al. 1992) indicated 
that the power spectrum was approximately scale-invariant (n ~ 1) on the 
large angular scales (£ < 30) to which COBE is sensitive, and this conclu- 
sion has been confirmed by the FIRS and TENERIFE experiments (Ganga 
et al. 1994; Hancock et al. 1994). Ground- and balloon-based experiments 
have produced numerous measurements of the fluctuations on degree scales 
{£ ~ 50 — 300), and although the situation is still far from clear, there is 
now some evidence that the power spectrum is larger in this ^-range than on 
COBE scales (e.g. Scott et al. 1995; Kogut & Hinshaw 1996). Recent results 
from the CAT experiment (Hancock et al. 1996; Scott et al. 1996) indicate 
that the power spectrum has fallen to lower values at £ ~ 600, which could 
be interpreted as there being a CDM-type "Doppler peak" around £ ~ 300. 

Our knowledge of the power spectrum for £ > 30 is still quite limited, 
since the small patches of sky surveyed at high resolution so far give a large 
sample variance. On the largest angular scales, however, the situation is 
much better, since COBE has mapped the entire sky. Indeed, the signal-to- 
noise ratio in the 4 year COBE data (Bennett et al. 1996) is so good that 
the error bars on the (^-estimates for £ < 10 are now entirely dominated by 
cosmic variance. This means that (apart from reducing possible systematic 
errors and correcting for residual foreground contamination), these are the 
best estimates that mankind will ever be able to make of the large-scale 
power. 

In view of this experimental progress, it is clearly worthwhile to estimate 
the power spectrum from the 4 year COBE data as accurately as possible. 
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This is the purpose of the present Letter. 

1.1 Power spectrum estimation with incomplete sky cover- 
age 

Pioneering work on the problem of power spectrum estimation from ex- 
perimental data (Peebles 1973; Hauser & Peebles 1973) has recently been 
extended and applied to the 4 year COBE data (Wright et al. 1996, here- 
after W96). When estimating power spectra, it is customary to place both 
vertical and horizontal error bars on the data points, as in e.g. Figure ^. 
The former represent the uncertainty due to noise and sample variance, and 
the latter reflect the fact that an estimate of Cg inadvertently also receives 
contributions from other multipole moments. In other words, the estimate 
of Cg is in fact a weighted average of a band of multipoles. The weights are 
referred to as the window function, and for the estimate of Ci to be a good 
one, we clearly want the window function to be centered on £ with an r.m.s. 
width A£ that is as small as possible. As is well-known, incomplete sky cov- 
erage destroys the orthogonality of the spherical harmonics, and makes it 
impossible to attain perfect spectral resolution, A£ = 0. For typical ground- 
and balloon-based experiments probing degree scales, the relative spectral 
blurring A£/£ tends to be of order unity, which makes it difficult to resolve 
details such as the number of Doppler peaks. A much better method is that 
used in W96, where the relative spectral resolution A£/£ is brought down 
to the order of 25% by using spherical harmonics. In Tegmark (1996a, here- 
after T96), a method was presented for reducing these horizontal error bars 
still further, down to their theoretical minimum, which for a 20° galactic 
cut was seen to be A£ ~ 1. 

1.2 The importance of high spectral resolution 

Just as high spectral resolution is crucial in for instance absorption line 
studies (since it prevents interesting features from getting smeared out), it 
is also important when measuring the CMB power spectrum. The reason is 
that we cannot a priori assume that the power spectrum will be a simple 
smooth function. Indeed, there are non-standard cosmologies, such as "small 
universe" models with nontrivial spatial topology, that predict power spectra 
which on average have ann~ 1 slope but contain bumps and wiggles that 
can only be resolved with a high spectral resolution (Stevens et al. 1993; de 
Oliveira-Costa & Smoot 1995). In addition, lowering the spectral resolution 
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degrades information: broad window functions make the estimates of nearby 
multipoles highly correlated, so that the resulting power spectrum plot will 
contain fewer independent data points than one may naively expect. 

The remainder of this Letter is organized as follows. In Section 2, the 
four year COBE data set is analyzed. In Section 3, the results are discussed 
and compared with those obtained with other methods such as the orthog- 
onalized spherical harmonic method of Gorski (1994, hereafter G94) and 
the signal-to-noise eigenmode methods of Bond (1995b, hereafter B95) and 
Bunn and Sugiyama (1995, hereafter BS95). 

2 RESULTS 
2.1 The method 

How to extract the angular power spectrum from a CMB map with maxi- 
mum spectral resolution is been described in detail in T96. We will apply 
this method to the 4 year data here exactly as it was applied to the 2 year 
COBE data in T96. Below we give merely a brief review of how the method 
works, referring the interested reader to T96 for technical details. 

A simple estimate of the multipole Cg is obtained by taking a linear 
combination of all pixels, squaring it, and subtracting off the expected noise 
contribution. The weights in the linear combination are conveniently plotted 
as a sky map in the same way that we plot the data, and we usually refer 
to the set of weights as a weight function. It is easy to show that the 
expectation value of such an estimate is a linear combination of all the true 
multipoles. The coefficients in this linear combination are called the window 
function, a function of I (see the examples in Figure ||). In other words, 
given any weight function, there is a corresponding window function. The 
window function turns out to be simply the square of the spherical harmonic 
coefficients of the weight function, summed over m. 

In the Hauser- Peebles method, the weight functions are chosen to be the 
spherical harmonics, but set equal to zero inside the galaxy cut and appro- 
priately rescaled. A number of other weight functions have been employed in 
CMB analyses, for instance the orthogonalized spherical harmonics of G94 
and the signal-to-noise eigenfunctions of B95 and BS95 (these functions were 
tailor made for the problem of efficient parameter estimation with likelihood 
analysis, not for power spectrum estimation). The method of T96 simply 
employs those weight functions that give the narrowest window functions 
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possible, and it is shown that these functions can be found by solving a cer- 
tain eigenvalue problem numerically. In the simplest version of the method 
(the version we use here — it has 7 = as in T96), these functions are com- 
pletely independent of any prior assumptions about the power spectrum, 
and are hence determined by the geometry of the galaxy cut alone. 

What are these weight functions like? A sample weight function is plot- 
ted in T96 together with the corresponding spherical harmonic, and it is 
seen that they look quite similar far away from the galactic cut. The main 
difference is that the optimal weight functions approach zero smoothly at 
the edge, whereas a truncated spherical harmonic does not. As is discussed 
in further detail in T96, this absence of sharp edges in the weight functions 
is a key feature of the method, since sharp edges cause "ringing" in Fourier 
space (in the multipole domain), which corresponds to an unnecessarily wide 
window function. 

To reduce error bars, a multipole is estimated as a weighted average of 
estimates of the above-mentioned simple type, just as in the Hauser-Peebles 
method where Ci is estimated by an average of the square amplitudes of the 
(2£ + 1) different spherical harmonic coefficients. 

2.2 The data 

The 53 and 90 GHz channels (A and B) of the COBE DMR 4 year data 
were combined into a single sky map by the standard minimum-variance 
weighting, pixel by pixel. We use the data set that was pixelized in galactic 
coordinates. After removing all pixels less than 20° away from the galactic 
plane, 4016 pixels remain. As has become standard, we make no attempts 
to subtract galactic contamination outside the cut. 

The resulting power spectrum is shown in Figure |J A brute force like- 
lihood analysis of the 4 year data set (Hinshaw et al. 1996) gives a best fit 
normalization of Q r ms,ps = 18.4 fxK for a simple n = 1 modelj^j correspond- 
ing to the heavy horizontal line in the figure. If this model were correct, 
we would expect approximately 68% of the data points to fall within the 
shaded 1 — a error region. As can be seen, the height of this region (the size 
of the vertical error bars) is dominated by cosmic variance for low I and by 

2 By this we mean an n = 1 model including only the Sachs- Wolfe effect, so that 
d oc Xjiil + 1). Note that the slow rise towards the first Doppler peak in an n = 1 
CDM model gives a best fit Sachs- Wolfe spectrum with n ~ f .15, whereas n = 1 models 
with spatial curvature or cosmological constant can give best fit Sachs- Wolfe spectra with 
n < 1. 
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Table 1: The COBE power spectrum 5T = [£(£ + l)Q/2vr] 1 / 2 in fxK. 

noise for large £. At the cost of increasing A£, the variance can of course 
be reduced further by grouping multipoles together in bands and averaging 
them with minimum- variance weighting, as shown in figures |3| and ||| and in 
Table 1. We have followed W96 and chosen the eight multipole bands 2, 3, 
4, 5-6, 7-9, 10-13, 14-19 and 20-30. 

For verification, 1000 mock COBE maps were generated for the model 
n = 1, Qrms,ps = 18.4 juK and piped through the data analysis software. As 
expected, the extracted multipoles were found to be unbiased estimates of 
the true multipoles, and the scatter was in agreement with the (analytic) 
error bars shown in the figures. 

2.3 The window functions 

The horizontal bars in Figure [l] are seen to be fairly independent of £, just 
as expected — as discussed in detail in T96, the angular size of the two 
sky patches surviving the galaxy cut is A9 ~ 1 radian, and we expect 
A£ ~ 1/A6 ~ 1. A typical window function is shown in Figure ^, and 
exhibits the following features that are common to all our window functions: 

• For even £, the window function vanishes for all odd multipoles, and 
vice versa. This happens because the galactic cut is symmetric about 
the galactic plane, and thus preserves the orthogonality between even 
and odd spherical harmonics, since they have opposite parity. 

• The window function is for all practical purposes zero for multipoles 
below £ — 2 and above £ + 2. 

• The central value is typically about three times as high as the two 
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sidelobes. 

The multipole Di = £(£ + l)Ci is thus typically estimated by something like 
0.6Dg + 0.2Di_2 + 0.2Z?£_|_2. This corresponds to all data points in Figure [l] 
being uncorrelated, with the exception that points separated by A£ = 2 have 
a correlation of order 55% and points separated by Al = 4 have a correlation 
of order 9%. Note that neighboring points are completely uncorrelated. 

The only exception to the above is the window function for the quadrupole, 
£ = 2: since it is required to vanish at £ = 0, it picks up a non-negligible 
contribution from £ = 6 instead. 

3 DISCUSSION 

We have computed the angular power spectrum of the 4 year COBE DMR 
data using the maximum resolution method of T96. The signal-to-noise ra- 
tio in the data is now so high that the error bars for £ < 10 are entirely 
dominated by cosmic variance. This means that, apart from future correc- 
tions due to better modeling of foregrounds and systematics, this is close to 
the best measurement of the low multipoles that mankind will ever be able 
to make, since cosmic variance could only be reduced by measurements in a 
different horizon volume. 

The power spectrum in Figure |l] is seen to be consistent with an n = 1, 
Qrms,ps = 18.4 fiK model. This model is close to the best-fit models found in 
the various two-parameter Bayesian likelihood analyses (e.g., Hinshaw et al. 
1996, Gorski et al. 1996, W96), which we can interpret as the best-fit straight 
line through the data points in Figure |] being close to the horizontal heavy 
line. As has frequently been pointed out (see e.g. White & Bunn 1995), 
Bayesean methods by their very nature can only make relative statements 
of merit about different models, and never address the question of whether 
the best-fit model itself is in fact inconsistent with the data. As an absurd 
example, the best fit straight line to a parabola on the interval [—1,1] is 
horizontal, even though this is a terrible fit to the data. It is thus quite 
reassuring that the power spectrum in Figure [l] not only has the right average 
normalization and slope, but that each and every one of the multipoles 
appear to be consistent with this standard best fit model.f] 



3 Kogut et al. (1996) find that correcting for galactic foreground emission increases the 
value of the CMB quadrupole, making it consistent with the best fit power spectrum. 
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3.1 Comparison with other results 

A number of other linear techniques for CMB analysis have recently been 
applied to the COBE data. Both the orthogonalized spherical harmonics 
method (G94), the Karhunen-Loeve (KL) signal-to-noise eigenmode method 
(B95, BS95), and the brute-force method (Tegmark & Bunn 1995, Hinshaw 
et al. 1996) were devised to solve a different problem than the one addressed 
here. If one is willing to parametrize the power spectrum by a small number 
of parameters, for instance a spectral index and an amplitude, then these 
methods provide an efficient way of estimating these parameters via a like- 
lihood analysis. Why cannot the basis functions of these methods be used 
to estimate the angular power spectrum Ce directly, as they are after all or- 
thogonal over the galaxy-cut sky? The answer is that these basis functions 
are orthogonal to each other, whereas in our context, we want them to be 
as orthogonal as possible not to each other but to the spherical harmonics. 
This is illustrated in Figure ||, which contrasts £ = 20 window functions of 
the optimal method and the generalized Hauser-Peebles method (W96, de 
Oliveira-Costa & Smoot 1995). We want the window function to be centered 
on I = 20, and be as narrow as possible, so the lower one is clearly prefer- 
able. The upper weight function is seen to couple strongly to many of the 
lower multipoles, and picks up a contribution from the quadrupole that is 
even greater than that from I = 20. This of course renders it inappropriate 
for estimating the power at I = 20. Analogous window functions can read- 
ily be computed for the orthogonalized spherical harmonics of G94 or the 
signal-to-noise eigenmodes of B95 and BS95. They are also broader than the 
optimal one in Figure ^ — the optimal weight functions of course give nar- 
rower window functions than other basis functions by definition, since they 
were defined as those functions that give the narrowest window functions 
possible. 

It should be emphasized that generating such window functions for the 
basis functions of G94, B95 and BS95 would be quite an unfair criticism of 
these methods, since this would be grading them with respect to a property 
that they were not designed to have. These authors have never claimed that 
their basis functions were optimal for multipole estimation, merely (and 
rightly so) that they were virtually optimal for parameter fitting with a 
likelihood analysis. 

The power spectrum of Figure [l] is also consistent with that extracted in 
W96 using the Hauser-Peebles method. At low t, the individual multipoles 
estimates agree well with each other. As t increases, the spectral resolution 
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A£ of the Hauser-Peebles method grows approximately linearly (see T96) 
whereas the resolution in Figure |l] is seen to more or less remain constraint. 
Thus the data points in W96 begin to differ from those in Figure [l] at larger 
£, since the former are no longer probing individual multipoles but weighted 
averages of a broad range. 

In summary, the 4 year COBE data has measured the power spectrum 
for i < 15 with an accuracy approaching the cosmic variance limit. As the 
next generation of CMB experiments extend this success to smaller angular 
scales, the CMB may turn out to be one of the most potent arbiters between 
cosmological models. 

The author wishes to thank Angelica de Oliveira-Costa, Krystof Gorski, 
Ned Wright and an anonymous referee for helpful comments on the manuscript 
This work was partially supported by European Union contract CHRX- 
CT93-0120 and Deutsche Forschungsgemeinschaft grant SFB-375. The COBE 
data sets were developed by the NASA Goddard Space Flight Center under 
the guidance of the COBE Science Working Group and were provided by 
the NSSDC. 
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Tegmark 1996 




Figure 1: The power spectrum observed by COBE. 
The observed multipoles Di = 1(1 + 1)C# are plotted with 1 — a error bars. 
The vertical error bars include both pixel noise and cosmic variance, and 
the horizontal bars show the width of the window functions used. If the 
true power spectrum is given by n = 1 and Q T ms,ps = 18.4^K (the heavy 
horizontal line), then the shaded region gives the 1 — a error bars and the 
dark-shaded region shows the contribution from cosmic variance. 
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Figure 2: Window functions before and after optimization. 
Two window functions for estimation of the multipole i = 20, m = 
are shown. The upper one is that of the spherical harmonic method 
(W96), which exhibits a strong leakage from lower multipoles such as the 
quadrupole. The lower one is the one resulting from the optimal method. 
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Figure 3: Power spectrum measurements. 
The observed multipoles Dg = £(£ + l)Cg are plotted for a selection of 
experiments. Both vertical and horizontal bars have the same meaning as in 
Figure |l[ The COBE data are those of Figure |], averaged over 8 multipole 
bands, and the rest are from the Saskatoon experiment (Netterfield et al. 
1996) and from the compilation of Scott et al. (1995). 
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Figure 4: Power spectrum models. 
The data points from Figure || with the narrowest window functions (COBE 
and Saskatoon) are compared with the predictions from four variants of the 
standard CDM model from Sugiyama (1995), all with n = 1 and = 0.05. 
From top to bottom at i = 200, they are a flat model with A = 0.7, a model 
with h = 0.3, the standard h = 0.5 model and a model with a reionization 
optical depth r ~ 2. 
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